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ABSTRACT 

We study the renormalization problem for the Hartree-Fock approximation of the 
0(A)— invariant 4 model in the symmetric phase and show how to systematically im- 
prove the corresponding diagrammatic resummation to achieve the correct renormaliza- 
tion properties of the effective field equations, including Renormalization-Group invari- 
ance with the one-loop beta function. These new Hartree-Fock dynamics is still of mean 
field type but includes memory effects which are generically nonlocal also in space. 



1 Introduction, summary and outlook 

Cosmology, particle and condensed matter physics have given, in recent years, a great 
impulse to the search of a deeper qualitative and quantitative understanding of out-of- 
equilibrium dynamics of quantum fields. In fact a treatment based on first principle of 
the late time and strongly coupled evolution of quantum systems would provide a better 
insight in an important class of phenomena such as the reheating of the universe after 
inflation or thermalization of the quark gluon plasma in the ultra-relativistic heavy-ion 
and hadron colliders (RHIC, LHC). 

The challenge in treating quantum field theories in non-equilibrium conditions is that, 
except for very short time, standard perturbation theory does not provide satisfactory 
results (this is true also at equilibrium with nonzero temperature). Therefore one has to 
look for nonperturbative approaches providing infinite partial resummations of Feynman 
diagrams |T| E] . The simplest of such schemes are the mean field approximations such as 
the leading-order large-N expansion j3J HJ El El and the Hartree, or Hartree-Fock (HF) 
variational method El El ED] • 
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These approximations have been extensively studied and their features are well known 
[TTJ [T3J EH EH UIH HHj : they do provide a backreaction term on the evolution of quantum 
fluctuations that stabilize dynamics at late time but, on the other hand, they fail to 
properly describe an important aspect of late time dynamics such as thermalization. More 
elaborate approaches going beyond mean field have been put forward by considering the 
2PI (or 2PPI) effective action ^3 ^] at two (or more) loop or at next-to-leading order 
in 1/N expansion [2*01 l2*Tj H5]. yielding indeed approximate numerical thermalization at 
strong coupling. 

Apart from the search for a better description of the late-time dynamics, there has 
recently been much progress in understanding more formal aspects of resummed approxi- 
mations such as their renormalization properties. This is a long standing problem [221 12H| 
and recently a systematic method has been developed [2HI23I2S1I2E] *° remove divergences 
in the ^-derivable approximations based on applying a BPHZ subtraction procedure to 
diagrams with resummed propagators. 

In this article we will consider the simple HF approximation of the O(N) 4 model 
as defined by variational principles, or equivalently by resummation of daisy (or bubble) 
diagrams. It is known [HI 1221123 that the usual renormalization of bare coupling and mass 
is not consistent, so that the simple Hartree approximation is not really renormalizable. 
We will show that this nonrenormalizability is due to the absence of leading logarithmically 
divergent contributions coming from diagrams which are not present in the standard HF 
resummation. We shall show that including also these contribution, plus suitably chosen 
finite parts, yields indeed a renormalized and Renormalization-Group invariant version of 
HF equations. We point out that these improved equations, although still of mean field 
type, are nonlocal in time (that is, there is memory) and generically nonlocal also in space 
space, unlike the original ones. We also verify that all these nonlocalities disappears in 
the N — > oo limit, which is well known to provide a renormalizable, RG invariant and 
local approximation to the out-of-equilibrium dynamics of the model. 

Section [21 is dedicated to generalities. In subsec. 12.11 the HF approach is introduced 
as a variational approximation with a Gaussian state functional in the Schrodinger pic- 
ture. The corresponding well known equations are rederived [see eqs. (12.5(1 and eq. ([2.7(1 ]. 
In subsec. 12.21 we review the CTP formulation of out-of-equilibrium problems and the 
equivalent definition of HF approximation as resummation of bubble diagrams. In sub- 
sec. 12.31 we establish some relations between the two approaches and define a generating 
functional T' [see eqs. ()2.20[) and ([2.21(1 ] whose proper renormalization would render finite 
the dynamical problem. 

In Section [21 we tackle the central problem of constructing renormalizable and RG 
invariant HF-like equations. In subsec. 13.11 we begin by studying, as a leading example, 
the static and homogeneous problem (i.e. the calculation of the effective potential) for 
the N = 1 case of a single scalar field. As already stated above, our analysis leads 
naturally to add contributions from diagrams absent in the HF approximation, in such 
a way to include all leading logarithmic divergences and construct a modified effective 
potential [see eqs. (|3.22() ] with the correct renormalization properties. In subsec. 13.21 the 
analysis of the HF approximation and its renormalizable improvement are generalized to 
the fully dynamical and inhomogeneous N > 1 case [see. eqs. ([3.48(1 ], Further remarks on 



2 



the freedom of choosing various initial conditions without spoiling renormalization and 
RG-invariance are made in subsec. 13.31 

There are several possible developments along the lines of this work. First of all it 
would be interesting to study numerically our modified HF field equations, to investigate 
how the space-time nonlocalities affect the time evolution as compared to the standard 
HF approximation, which is known to fail even qualitatively at late times. Secondly, with 
the proper changes in the renormalization scheme, the derivation could be extended to 
the case of broken O(N) symmetry, where the standard HF approximation faces other 
infrared difficulties connected with the Goldstone theorem Since our improved HF 
resummation adds also cutoff-independent contributions to the conventional one, it is 
conceivable that infrared properties of the model are affected too. Another challenging 
task is the extension of our approach to the full two-loop 2PI effective action, through the 
inclusion of the nonlocal sunset diagram which is absent by definition in any mean field 
approximation. In fact one should expect that, even if such an inclusion allows to recover 
renormalizability as compared to the conventional HF approximation, the two-loop 2PI 
self-consistent equations still lack RG invariance with the two-loop beta function, since 
the 2PI effective action does not contain all diagrams which contribute to the next-to- 
leading ultraviolet divergences. 

2 Generalities 

A non-equilibrium approach to quantum field theory is needed every time we have to deal 
with an initial value problem. In fact, the usual in-out formalism of QFT provides the 
means to calculate scattering expectation values 

(in\0(t 1 ,t 2 ,...)\out) 

where \in) and \out) represent particle states prepared at a distant past and future, 
respectively, while O is some observable depending on intermediate times ti,t 2 , . . . . On 
the other hand, in order to study the real time dynamics from a given initial condition 
we need to know amplitudes like 

(*(t0)||O(tl,t2,...)l*(*0)) 

where is a generic state prepared at a initial time to < t\, t 2 , ■ ■ ■ ■ Out-of-equilibrium 
QFT provides the general setup for the calculation of such matrix elements, as well as 
more general expectation values in statistical mixtures of pure states like \^). 

In this section we briefly review some generalities about non-equilibrium QFT and 
point out some properties that will be useful to derive the results of the next section. 

2.1 Variational approach and HF approximation 

A simple and intuitive approach to initial value problem is to treat QFT as ordinary 
Quantum Mechanics of extended systems. We restrict here to the case of interest for this 
paper: scalar field theory in 3 + 1 dimensions with quartic interaction and unbroken O(N) 
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symmetry. The field variables are <fi(x) where x = (xi,X2,x^) are space coordinates and 
i = 1, . . . , N is the O(N) index. Classical dynamics is defined by an action, functional of 
trajectories fi{x), with x = (x, t), in the form 

S[(p] = J d 4 x ^d^i{x)d^^i{x) - ^-^(x)^(a;) - ^[^(x)^(x)] 2 | (2.1) 

where m 2 and A are the squared mass and the coupling constant respectively. In terms of 
(Pi(x) and its canonical conjugated momentum n^x) = (pi(x) the Hamiltonian reads 



H[cp, tt] = J d 3 x 



\i(x)iXi(x) + -V<fi(x) ■ V^(a;) + V(<fi(x)) 



where 

Canonical quantization proceeds by imposing the standard commutation rules 

[(fj{x), itk{y)) = i5 jk 5 3 (x - y) 

In the "position" representation states are functionals of ipi(x) and the momentum oper- 
ator reads 

7Ti(x) = -I 

5ipi(x) 

It is then straightforward to write down the Hamiltonian operator and the corresponding 
Schrodinger equation. 

In this framework, the standard variational methods of Quantum Mechanics suggest an 
easy nonperturbative approximation to dynamics (see e.g. ^01 EDI HH] ) • The variational 
method follows from the observation that Schrodinger equation is obtained by minimizing 
the Dirac action 

Sb[*]= / dt (*(t)\(id t -H)\*(t)) 

with respect to trajectories of vector states |^(t)) in the Hilbert space. Therefore minimiz- 
ing Sb over trajectories in a chosen variational family of states |^ V ar) gives an approximate 
solution of Schrodinger equation. In QFT, where scalar products are written as functional 
integrals, calculability strongly reduces the allowed family of states. Choosing states rep- 
resented by Gaussian wavefunctionals, for which calculability is manifest, yields the HF 
approximation. Let us in fact consider the wavefunctional 

*[<p](t) = M exp {i (p(t)\p - <f>{t)) -{<p- <j>(t)\ [\g~\t) + iS{t)] \<p - (f>(t))} (2.2) 

where the variational parameters are the background field <f>(x), the background momen- 
tum p(x), the real symmetric positive kernel Qij(x,y;t) and the real symmetric kernel 
Sij(x, y; t) with the short-hand notation 



(a(t)\M(t)\b(t)) = / d 6 xd 6 y aj (x,t)M jk (x,y;t) b k (y,t) (2.3) 
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By evaluating Sd we obtain the variational action for these parameters 



r HF [^p,g,s] = / dt((p(t)\<j>(t)-p(t))-v[g(t),<j>(t)] 



+ Tr 



g{t)s(t) - 2S{t)g{t)s(t) - \gity 1 



(2.4) 



where traces are taken over all indices and space variables, and 

V[g,<t>] = ^(<f>\-A + ml\<P) + ^Tr[(-A + ml)g]+ j d 3 x( 



4! 



+ ^(p i (x)(pi(x)g jj (x,x) + ^4>i{x)4)j{x)gij{x,x) 

+ ^gii(x, X)g jj (x, x) + -^pijix, X)gij(x, x) 

Notice that the HF Dirac action is a first order action that involves only the first time- 
derivative of fields and treats coordinates and momenta as independent variables. The 
equations of motion are obtained by variation of eq. (|2.4j) as 

<Pi = ~ { [—A + ml + |A 4>k4>k] $ij + f A T ijkm g km (x, x)}(f)j 

g iS = 2[g s + s g}ij 

S l3 {x,y) = [\g- 2 -2S%{x,y) 

- \ [(-A + ml )5ij + ±A r ijkm ((j) k (j) m + g krn (x, x)] 5 3 (x - y) 

where 

7~ij k m ^{$ij$km ^ik^jm $im$jk) 

and matrix multiplication, over both discrete and continuous indices, is understood. We 
have tacitly replaced in the last equations, and in general in the quantum Hamiltonian, 
the constant parameters A, m 2 with bare (cut-off dependent) parameters Ao, m\. In fact 
the theory should be thought as regularized with an UV cut-off A and then renormalized 
to remove the divergent dependence on A. We have not included here any field renormal- 
ization because, as is well known, this is absent in any mean-field type approach such as 
the HF approximation. 

To conclude this subsection we introduce an equivalent formulation of the eqs. (j2.5|) 
in terms of mode functions Uk a {k is the wavevector and a the O(N) polarization), that 
will be useful later on. For simplicity let us suppose that the initial (t = 0) kernels 
are transitionally invariant (while the background field and momentum may be point 
dependent). We can then write 



Gij(x,y; 0) 
Sij(x,y;0) 



d " k a tJ (k) e ikx 



(27r) 3 " y 
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Next, we introduce the t = mode functions by 



u kai (x,0) = \g(k) 



1/2] ikx 



U k a,i{x,0) 



^(fc)- 1 + 2S(k) u kaJ {x,0) 



(2.6) 



and let them evolve according to the equations of motion 

d 3 p 



(□ + ml) Sij + \\ t, 



t]mn 



[X 



ml* rrn 



X + 



(2tt)3 p 

Then one can show that the last two equations in eqs. ()2.5j) are equivalent to 

d 3 k 



Uka,j{x) = 

(2.7) 



Qij(x,y;t) = J u kaji (x,t)u kaJ (y,t) 



{21 



u ka ,i(x,t) = J d 3 y<^--[Q 1 } ij (x,y;t) + 2Sij(x,y;t)>u k a,j(y,t) 

Through the mode functions we can introduce the symmetric correlation 

d 3 k 



G i:j (x,y) = Gji(y,x) = Re J 
whose equal time value reproduce the kernel Q, 



(27T) 



(y) 



Gij(x,y;t) 



V\ m'\x =y =t 

Then we can reformulate the dynamics in terms of and G in a manifestly covariant way 



as 



{ [□ + ml + |A (f) k (x)(j) k (x)] + |A r ijkm G km (x, x)}(j) j (x) = 
{(□ + ml) 8^ + \\QTij km [<j>k{x)(j) m {x) + G km (x,x)] } Gij(x,y) = 

2.2 CTP formalism and resummations 



(2.9) 



The general approach to non equilibrium problems in QFT was developed by Keldish 
and Schwinger. It is known as closed time path (CTP) formalism and allows to use 
standard path integral functional methods (see [T71 [TBI E2 EH])- Basically it is obtained 
by introducing path integrals on a time path going from t = to t = +oo and back. Field 
integration variables are then doubled and subdivided into (+)— components, for the path 
integral forward in time, and (—)— components for the backward piece. Given an initial 
state defined by the functional ^[<f\ (in our case it will be a Gaussian state) one writes 
down the functional integral 



iW\j+,j- 



iS{<p+]-iS{ip-]+i{j+\<p + )~i(j-\ip-) 



(2.10) 
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where S is the classical action (|2.1|) and we have used the notation ()2.3|) for the currents 
terms. Integration is on trajectories from t = to t = +00 (with the condition ip + = ip_ 
at t = +00) and tp± in the wave functional is the t = section of <p±. By construction 
W[j + , j-] is the generating functional of connected Green functions 



Vn) 



> 



n+m rn+m 



5j + {xx) . . . 6j + (x n )6j-(y!) . . . 83 -(v* 



j+=o 

j-=0 



= -i {*\T{<p(y x ) . . . ^( 2/m )}T{^(a; 1 ) . . . <p(x n )} |^) conn 

where T and T define time ordered and inverse ordered products, respectively, and in- 
ternal indices have been omitted for ease of notation. The effective action r 1PI [</> + , </>__], 
which is the generator of 1PI vertex functions, is the Legendre transform of W[j + ,j_] 
from the currents j± to the fields 0±. The equation of motion for the background field 
<f>(x) = ^{x) \^>) then reads 



5T 



ipi 



5(f) + (x) 







(2.11) 



Notice that, in the case of the Gaussian wavefunctional eq. (|2.2|) . Tipi depends paramet- 
rically only on the t = kernels Q{0) and S(0), while the t — background field 0(0) and 
p(0) enter, instead, as initial conditions for the equation eq. (|2.11|) . 

The perturbative diagrammatic expansion in the CTP formalism proceeds as in vac- 
uum QFT in terms of free propagators 



G ( LJx,y) = G° F (x,y) 



' + 

G°_ 



ix,y) 



G° P (x,y) = 
G°_Jy,x) 



-i(^\T<p{x)ip{y)\^) 
-i (*| T<f{x)<f{y) I*) 
-i(V\<p(y)<p(x) |^) 



conn 



I A=0 



conn 



A=0 
conn IA=0 



and vertices (for clarity, we write here explicitly the internal indices) 




i X 



T ijkm 




T ijkr. 



(2.12) 



It can be shown that, in generic out-of-equilibrium contexts or even at equilibrium with 
nonzero temperature, plain perturbation theory is of very little help and many resum- 
mation methods have been developed to go beyond it. A very successful instrument in 
this sense is the 2PI effective action (see [13 CHI)- It is defined as the double Legendre 
transform of the W generating functional with respect to the usual current one-point j± 
and to the two-points current K a p{x, y) coupled through the term 



d 4 x / d 4 y K a p(x,y)Lp a (x)tp ( 3(x) 
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where a, (3 = ±. r 2 pi is a functional of the classical fields a and of the propagators G a p. 
It yields two equations of motions 



5Y 



2PI 



<5r 2PI 



S(j) Q (x) * ' 5G a p(x,y) 



(2.13) 



Here the notation I* indicates that, by their physical meaning, the (±)-component fields 
and propagators have to satisfy, on the solutions of motion, the following relations 

4>-{x) = 4>-{x) = 0(x) 
G F (x, y) = G+-(y, x)0(x o - y ) + G + _(x, y)9(y - x ) 
Gp(x, y) = G+_(x, y)9(x - y ) + G + ^(y, x)6(y - x ) 

Hence the system fl2.13j) reduces to two coupled equations for <fi and G+_ only. Moreover, 
any initial Gaussian state may be absorbed in the t — term for the j a and K a p currents, 
so that, by the double Legendre transform the initial Gaussian state disappears from 
the effective action, but fixes the initial conditions on <p and G + _. The role of 0(0) and 
p(0) = 0(0) is immediate, while for the kernels we have 

G+-(x, y)\ Xo=yo =o = -% (*| <p(y)<p{x) |^) - i (f>(y, 0)<f>{x, 0) = Q(x, y; 0) 

d 

—G + _(x,y)\ xo=yo=0 = -i(y\ir(y)(f)(x) |*) - ip(y, 0)<f)(x, 0) 

= 2i [gS}(x,y;0) + ±5 3 (x-y) 

Given r 2 pi at a certain perturbative loop order, if we solve the second equation in ()2.13|) 
for a generic <fi and substitute the result G[(p] into the first one we obtain the background 
equations of motion corresponding to a resummed diagrammatic approximation of the 
1PI effective action Tipi. 

For a scalar theory, the T 2 pi has the general form 

r 2P i [0, G) = S[</>] + ^Tr [logG] + ^Tr [G^G] + T 2 [0, G] 

Here S is the complete classical action of the double time path (i.e. S = S + — S-). 
Traces are taken over all indices i, a and x. Gq 1 is the second derivative of the action 
in a background, T 2 is the sum of all vacuum 2PI diagrams with G propagators and 
vertices defined by the classical action in a background. To two loops level the diagrams 
contributing to the T 2 are the "8" and "sunset" diagrams 




(2.14) 

The HF approximation corresponds to consider only the first contribution to T 2 . 

r 2 = |A [G 2 F {x,x)-G 2 p (x,x)] 
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In fact using this form of T in the field equations (j2.13j) . setting G(x,y) = y) + 

G- + (y,x)} and observing that the antisymmetric combination decouples, one obtains 
exactly the HF equations (j2.9|) . 

Notice that "8" is the only 2PI diagram made of "product" of loops corresponding 
to a mean field contribution to the mass. In the 1PI framework this corresponds to a 
resummation of all vacuum 1PI diagrams with daisy and superdaisy topologies of the 
form 




2.3 Physical representation and reparametrization of HF effec- 
tive action 

In this subsection we derive useful explicit relations connecting the diagrammatic defini- 
tion of Thf to the variational approach of subsection 12.11 To this purpose it is convenient 
to introduce a different representation of the CTP formalism, known as the physical rep- 
resentation (see [US!)- We introduce the field redefinitions (omitting again the internal 
indices to simplify notation) 

0A = 0+-0- , c =^(0+ + 0_) 



and write the 1PI effective action as a functional of these new fields, Tipi 
calculating vertex functions, one then finds 



By 



8 n T 



(xi) . . . 5(j) c (x n ) 



and 



8 



n+mi 



5(f) c (xi) . . . 5cj) c (x n )5(j>A(yi) ■ ■ ■ S<j) A (y ri 







(2.16) 



(2.17) 



0A=O 



if the time component of anyone of the x's is larger than the time component of anyone 
of the y's. Let us also remark that, by the definition of CTP generating functional, all 
time coordinates in the vertex functions are supposed to be positive so we can set, as well, 
these functions to be zero for any negative time. Using eq. ()2.16|) together with eq. ()2.11|) 
one can write the equations of motions in the form 



5T 



(2.18) 



i>A=0, <f>c=<t> 
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Notice that 2n-legs vertex functions with one <p A leg and In — 1 C legs are the only ones 
contributing to these equations of motion. Then eq. (j2.17J) guarantees that all the terms 
nonlocal in time in eq. (j2.18j) do satisfy causality. 

Perturbative calculations by diagrammatic expansion in the physical representation 
are based on the free propagators 



G° cA (x,y) 
G° Ac (x,y) 



G° A (x, y) = -i6{y Q - x ) (¥| [0(x), cj>(y)] |*) | A =o 
G° R (x, y) = -i9(x - y ) [<f>(x), <t>(y)] |tf > | A =o 
-2i (*| {<f>(y), <f>(x)} \^) conn | A=0 , G° cc (x, y) = 



The free retarded and advanced Green functions G° A and G° R do not depend on the initial 
state and are translational invariant. The correlation function G 
depend on The vertices are 



iG AA , instead, does 




i A r, 



ijkm 



^0 T ijkr. 



(2.19) 



where solid lines represent </> c legs while dotted lines represent (f) A legs. 

As stated above, the HF approximation consists in the resummation of diagrams with 
daisy and superdaisy topologies. Therefore the functional differentiations which produce 
the vertex functions act only on propagators; hence the diagrams contributing to a 2n-legs 
vertex function (both in ± and physical representations) must have the form 




that is, the In legs are pairwise connected to to n vertices. This means that the HF 
effective action can be written in the form 

r H F[<M,0 c ] = - (<MI □ 1 C ) - FlCxiv] 

where T is a functional of the following composite matrix fields 

iij[x) = C) i(x)0 c j(x) 
Xij(x) = c ,i(x)0 A j(x) 

r)ij(x) = (f) Aji {x)(j> A> j(x) 
We now introduce, for ease of notation, the new object 

1 bT 



2 5xtj(x) 
10 



x=v=o 



which is a functional of £ij(x) = 4> c ^(x)(f) C} j(x) = <pi(x)(pj(x) only. Then the equation of 
motion in the HF approximation takes the general form 

(□^• + 2^)^ = (2.20) 

By comparison with eq. ()2.5j) and eq. (|2.8|) we obtain the following expression for T' 



(27T) 



3 V'pa, i\<£ ) ^pa,j\^) 



(2.21) 



The dependence of the mode functions u pa on £ is fixed by solving eq. (|2.7|) for a generic 
field. Notice that the identification of eq. ()2.21|) is somewhat arbitrary since we can add 
to T' a generic functional term A[£\ icky (x)£, km (x) with A icky = —A skim without changing 
the contribution of T' to the effective action and to the equations of motion (J2.2U)) . 

The role of the functional T' in the generation of vertex functions can be easily estab- 
lished. Let us consider the 2n-points vertex function with one 0a field leg and m(= 2n — 1) 
4> c field legs (as we said those we are interested in for they contribute to the equations of 
motion), a simple calculation leads to 



Vm) 



S 2n T 



5<pA i(x)S(p c h (y{) ...5(j) c jm (y m ) 
□ 5 A (x- y jl )5 1>n + 

n-1 

Y[^(py 2k - PV2k+i) 



[n 



k=l 



</>c=0 



— S 4 (x-py 1 ) 



(2.22) 



Where S m is the set of all permutations on {1 . . . m} and we used the notation pyi = y p ^, 
P3i = 3 P (i) and /<,,;, i//.) = C,,;,.^ (/'//>)• 

In what follows we are going to study the problem of renormalizability (and RG- 
invariance) of the effective field equations in the HF approximation. By the results of 
this section we need to consider only the vertex functions contributing to the equations 
of motion and hence only the functional J 7 '. 



3 Renormalizability and RG invariance 

We are concerned here about the possibility of constructing a set of HF equations in 
which all the divergent dependence from the cut-off has been removed by a suitable 
renormalization procedure and, at the same time, no dependence on renormalization 
scale has been introduced. In brief we are going to study the problem of renormalizability 
and RG invariance of HF equations. This will lead us to some results that go beyond 
the strict mean field approximation and include resummations coming from 2PI effective 
action loop expansion. 
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In this section we will deal with the general case of an inhomogeneous (in space) dy- 
namical problem within the O(N) symmetric scalar field theory with unbroken symmetry. 
However the main aspects of the problem can be more easily pointed out by considering 
a simpler case. In the first subsection, therefore, we restrict to the static problem, that is 
the calculation of the effective potential, for the N = 1 theory. 



3.1 A simple instructive example: effective potential for N=l 

We want to study the renormalization problem of the static and translationally invariant 
version of eqs. (|2.5|) when N — 1. This situation follows, in the variational approach of 
sec. 12.11 by considering a field <ft constant throughout space-time, with vanishing momen- 
tum p, and a constant translationally invariant Q kernel 

d 3 k 



G{x ~ y)= S m e ~ lk{x ~ v) 



with S = 0. By a simple substitution in eq. (J2.5|) we obtain 

1 1 2 i 2 , ,2 , f d 3 p ~ 



4£ 2 (fc) U 2 2 7 |p|<A (2vr)3 





2 , ^0 , 2 , ^0 f d 3 p 



m ^ + -ji^ + -jn / -^-g( P ) 

6^ 2 7 |p|<A (2vr)3 yw 



(3.1) 



where a sharp cut-off A has been introduced as regularization. 

Equivalently we can consider the equations of motion as defined by the 1PI effective 
action and reduce to the static case by evaluating r 1PI on a space-time constant 0. So 
letting 

we can define the effective potential by Kff(0) = V(<p 2 ), with V(£) the primitive of V'(£) 
vanishing at £ = 0. From eq. ()2.20j) the equation for the background field <fi is 

K' ff (0) = 2H0 2 )0 = o 

Comparison with the second equation in 1)3.1)) leads immediately to 

where the implicit dependence of Q on £ = 2 is obtained by solving the first eq. ()3.1|) 
with a generic 0. We recall also that V e s(4>) is the generator of vertex functions with all 
incoming momenta set to zero (for constant homogeneous fields the distinction between 
CTP and vacuum standard formalism is immaterial). 
As usual, one introduces the ansatz 

~ 2Vk 2 + M 2 
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that allows to cast the first of eqs. (j3.1|) in the form of a mass gap equation 



2 , A ^ , A f drp 



^ ( f ) = ^ + f£ + f y^^__ p.3) 

where the integration has been written in four Euclidean dimensions (p 2 = pi + p 2 ). By 
comparing eq (|3.3jl with eq. f)3.2j) we can rewrite the latter as 

f.(0 = ^_£t (3.4) 

so that, by combining eq. ()3.3|) and eq. (|3.4jh we have 

^(0 = <+^+h [ pL. 5 1 (3.5) 

i ; 2 12* 4 7 p2<A2 (2tt)V + 2V'(£) + §A £ 

This is a self-consistent definition of the effective potential and is the central equation 
that we are going to use in this subsection. 

Let us apply the standard renormalization procedure. We introduce two physical 
parameters, the renormalized squared mass m 2 and the coupling constant A, which are 
identified, respectively, as the second and fourth derivatives of the effective potential at 
zero field, so that 



d 2 V r 



cff 



0=0 



2V'(0)=m 2 , =12V"(0) = \ (3.6) 

0=0 



These are the two renormalization conditions that, once inverted, determine the depen- 
dence of the bare parameters Ao(A, m, A) and m^(X, m, A) on the physical ones and the 
cut-off. The assumed existence and positivity of 2V'(0) provides the definition of unbroken 
symmetry, which is our choice here. 

Using eq. (|3.5|) . we can calculate the explicit HF form of eq. ()3.(i|) 



m 2 = m 2 + ^>(m 2 ,A) , A= - 2A (3.7) 

z l+2Ao^(^ 2 ,A) 

where, to shorten notation, we have introduced 

I%\m 2 ,A)=[ ^ V - 1 (3.8) 
J p2<A 2 (2vr) 4 (p 2 + m 2 ) n 

and the suffix refers to the Euclidean form of the integral. As A — ► oo one sees that 
Jg is quadratically divergent in A, while 1^ is logarithmically divergent and all other 
1^ are finite. 

We now point out the problems that occur with this standard procedure. 

First of all we can see that there is a pathological behaviour of the effective quartic 
coupling A as a function of the bare parameter Ao at fixed cut-off A. In fact A has the 
correct 1— loop A 2 , term dictated by perturbation theory, but certainly fails at higher 
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orders, since is exhibits an unphysical behaviour, growing to a maximum value at Ao = 
\max^ then decreasing to zero and to even more unphysical negative values. This implies 
the breakdown of the HF approximation for values of Ao greater than \™ ax in a theory 
at fixed cut-off. Another way of looking at this problem is to invert eqs. (|3.7|) to obtain 
Trip and Ao as functions of A parametrized by m 2 and A. A full trajectory in A of TOq and 
Ao corresponds to a single renormalized theory describing the dynamics at momentum 
scales much smaller than A. Clearly we must restrict our attention to the monotonically 
increasing branch through the origin at Ao = A = 0. But these trajectories, for any 
positive value of A, exhibit a "Landau obstruction" at the value of A which corresponds 
to \™ ax . This obstruction is even more troublesome than the Landau pole present in the 
standard 1— Loop-Renormalization-Group improved relation, which by a suitable choice 
of finite parts can be written as 

A ° (A ' A/m)| — = l-lA/f(^,A) <3 - 9> 

in terms of the whole Ig (m 2 ,A) rather than just its leading divergence ^2 log(A/m). 
In fact, the obstruction spoils even the one-to-one correspondence between bare and 
renormalized parameters (at fixed cut-off) which must hold true in general and holds 
true also in the 1LRG improved relation ()3.9|) . 

Secondly, and even more seriously, one may verify that the above renormalization 
procedure does not remove all the logarithmic UV-cutoff dependence from the potential. 
That is, V e g(<f)) does not parametrically depend solely on m 2 and A, but also on Ao and 
therefore on log A. In a QFT with a Landau pole or obstruction, where the UV cutoff 
cannot be completely removed, we should accept at most an inverse power dependence of 
physical quantities on A. 

To show this, consider the finite part of 1$, that is 

J(y, A) = I§\y, A) - 4V 2 , A) - 4 2) (m 2 , A)(m 2 - y) 
for any real y. By the use of J eq. (J3.5|) may be cast in the form 

V'(0 = \m 2 + ^A£ + -L(A + 2 A ) J (2 f '(£) + ±A £ - m 2 , A) 
This can be written more compactly 

*V'(t) = \gJ(2AV , {t) + lgt) (3.10) 

with the definitions 

S = !(A + 2A ) , AV'(O = V'(O-V'(0)-V"(0)£ (3.11) 

It is quite clear now that the higher order derivatives (or vertex functions with vanishing 
incoming momenta) of V e s, which are all contained in AV f , depend explicitly on logarithms 
of the cut-off through the effective coupling g. 
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This problem was dealt with in [SJ |2H] by a direct renormalization of HF gap equation, 
by setting 

M 2 = m 2 + ±A £ + |A 4 1} (M 2 , A) 



m 2 + \\£ + \\ J(M 2 ,A) 



so that the bare-to-renormalized relations are found to be 

Ao= n 4- - , m 2 = m 2 -iAo4 1} (m 2 ,A) (3.12) 
1 - ^Xr E '{m 2 ,A) 

This procedure makes M 2 (£) and all higher derivatives of V(£) cutoff independent by 
definition, but fails in two respects. First of all the 1— loop beta function read out from 
eq. (|3.12|) is three time smaller than the correct one [see eq. ([3.9)1 ] . which implies a mis- 
match already to first order in perturbation theory. Moreover the parametrizations (|3.12|) 
do not eliminate completely the cutoff logarithms, which remain in the quartic term of 
Ks(0) since the background field equation reads 

VM) = [|M 2 (0 2 ) - |A O 2 ]0 = (3.13) 

One could adjust the mismatch with the beta function by redefining the bare quantities 
appearing in the HF gap equation. If we make the ad hoc replacements 

Ac 



An ► Ar 



l + Ao4 2) (m 2 ,A) (3.14) 



2 k ~2 _ 2 l\ rWf m 2 



m n — > m n = m 



iA Q /W(m 2 ,A) 



then the correct 1LRG relation ()3.9j) is recovered. As for the problem with the quartic 
term in V e g((f)), a reasonable argument was put forward in [B] to justify the substitution 
An ->• A in eq. (jSH). 

We will see that this substitution, as well as the redefinition in eq. ([3. 14)1 . can be fully 
justified in the spirit of the 1LRG improvement of the HF approximation by an explicit 
diagrammatic analysis. 

The problems just outlined, have, in fact, a simple interpretation if we analyze the 
diagrams exactly resummed by the HF approximation. First, if we consider the renor- 
malized mass, given by the 2 legs vertex function with no momenta V"'(0), we see, by the 
second equation in (|3.7|) that it is dressed by daisy and superdaisy tadpole diagrams, as 





(3.15) 

From now on we are going to work with this dressed mass and we will draw a simple solid 
line the corresponding dressed propagator. Let us consider the quartic coupling, the first 
equation in ()3.7|) tells us it is obtained by "chain" diagrams (with dressed propagators). 



X 
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This is not the full diagrammatic expansion that corresponds to the 1LRG relation be- 
tween A and Ao in eq. (|3.9j) . The latter in fact includes all local contributions with Leading 
Logs (LL) of the cutoff from all diagrams at every given loop order. So eq. (j3.16J) misses 
the LL local parts of graphs like 




sym 




sym + 



(3.17) 



This is the cause of the problems outlined above with the form eq. (|3.7|) of the renormalized 
coupling A. 

Moreover, we can examine in the same spirit the diagrams contributing to higher order 
vertex functions (higher order derivatives of V in zero field). From eq. ()3.5|) we see that 
they are convergent loop diagrams (loop with more than 2 propagators) with dressed 
propagators 







(3.18) 



and vertices with the effective quartic coupling g = (2Aq + A)/3 



X 



(3.19) 



The reason of the presence of g and not of the renormalized coupling A in higher order 
vertices is that the local LL contributions of diagrams like 




(3.20) 



are not included in the resummation defined by the standard HF approximation as de- 
scribed by eq. ()3. !())) ■ Including such LL contributions would make the HF approximation 
renormalizable in terms of the reparametrizations of eq. ()3.7j) . that is with the patholog- 
ical relation between A and Aq. This is because the LL contributions of diagrams like 



(3.21) 



are still missing. Their proper inclusion would allow to obtain everywhere the renormal- 
ized 1LRG improved coupling constant. 
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By these observations we are led to a natural way to modify the HF approximation in 
order to obtain a renormalized effective potential and the form (|3.9|) for the coupling con- 
stant, that will ensure RG invariance in the general case, as we will see. The prescription 
is to include all missing contributions of the type (|3.17|) to 12 V"(0) in order to obtain 
the 1LRG relation between A and Ao as in eq. (|3.9|) . and to include also the contributions 
of graphs of the form (|3.2U|) and (|3.21j) in order to obtain A instead of g in higher order 
vertex functions. 

In conclusion we define a modified version of HF approximation by setting 

2V'(0) = m 2 = m 2 + \\ I^\m\K) 

12V"(0) = A = -— ^ (3.22) 

l + |Ao4 2 V,A) 

AV'(0 = JA J{2/\V\i) + \\() 
This defines to the renormalized potential V(0) and leads to the following equations 

M 2 = \\i + m 2 + \\J(M 2 ) 
O=(|M 2 -|Ae)0 

these are finite and exactly in the form found in [H]. This gives a simple diagrammatic 
explanation of the substitution Aq — > A made in jH] and of the statement (|3. 14)1 . 



3.2 The general case 

We are now in a position to deal with the general renormalization problem of the time- 
dependent and inhomogeneous HF approximation for O(N)— symmetric theory, as de- 
scribed by eqs. fjOjl . (|2~7jl . 

Our first step is to derive a self-consistent functional equation for the effective action 
analogous to the the one for the effective potential in eq. (|3.5jl . Let us introduce the 
squared mass parameter m 2 . We choose m 2 to be the equilibrium mass (i.e. the physical 
mass of the unbroken symmetry phase), that is to say the solution of the mass-gap 
eq. f!3.3|) . generalized to iV > 1, in zero field. 

m 2 = m 2 + i(AT + 2)Ao4 1) (m 2 , A) (3.23) 
Let us then define the free mode functions u ka as solutions of 

(□ + m 2 )4°i(x) = (3.24) 
with the same initial conditions of the exact mode functions, that is 

u < £ ) a (x,0) = u ka (x,0) , v$l(x,0)=u ka (x,0) (3.25) 
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We need also the free retarded and advanced Green functions 



(?) ^ _ „\ - f d p 6 Ji e -»P(*-») 



(27r) 4 p 2 — m 2 + iepo 



/?(0) c_ _ „N _ / _j*J! 5 V p -i P (x-y) 

A,ij \ y> J ( 2vr )4 J9 2_ m 2_, epo 

which allow to cast eq. (|2.7|) in the form of an integral equation (this is a quite standard 
way to proceed in the renormalization of out-of-equilibrium problems, see i.e. |34| ) 

«fca,iO) = uf a>i {x) + J d 4 yGf in (x - y) V nj (y) u kaJ (y) (3.26) 

where 

Vij(x) = |A ^j(x) + 2 r [e] y (ar) - m 2 % (3.27) 

plays evidently the role of mean field for the mode functions. In these equations, and 
everywhere else from now on, all field (0, mode functions, J 7 ', V, etc..) are to be thought 
as defined only for positive times (initial conditions are at the limit point t = + ) and 
all time integrations are restricted to positive values, as appropriate in an initial value 
problem. 

In compact notation eqs. (j3.26j) . ()3.27|) can be written as 

u ka = uf a + G^Vu ka , V = |A £ + 2r-m 2 (3.28) 

where sums over internal O(N) indices and integrations over spacetime coordinates are 
understood. In particular, notice that V, £ and act as multiplication operators over 
spacetime. 

The solution of eq. ()3.28|) formally reads 

= [1 - GSVr 1 ^ 

so that the cutoffed correlation 

Gij(x,y) = Re / — u pai (x) u pa j(y) 
J\p\<a ( 27r ) 



can be written 
where 



G = G[V] = [1 - G^V}' 1 G (0) [1 - V T Gf}~ 1 (3.29) 

G^ix -y)=Re J |^ uf a ,{x) 

is the free correlation function, which is entirely fixed by the initial conditions. In conclu- 
sion, by substituting, the form (j3.29j) in eq. (|2.21j) . we obtain the sought self-consistent 
equation for T' as 

( (i) A (3 ' 30) 



where we have renamed the correlation at coincident points as 



I[VUx) = G[VUx,x) 

Now, following closely the previous subsection, we may use these expressions, together 
with the general formula eq. (|2.22|) . to evaluate the renormalization conditions. To begin 
with, we consider the HF vacuum as initial state for the quantum fluctuations. That is to 
say that we start from equilibrium initial conditions for the mode functions by choosing as 
t — kernels the solutions of the static problem of the previous subsection (generalization 
to iV > 1 case is straightforward), that is 

^(fc) = ^y, u(k) = Vfc 2 + m 2 , 4-(*) = ( 3 - 31 ) 

Let us stress that, in spite of this choice, we are still considering an out-of-equilibrium 
problem since we allow for generic initial conditions for the background field. We will 
discuss the use of different initial kernels later on. 

By the standard prescription we have to compute, using eq. (|2.22|) and eq. (|3.30jl . 
(the Fourier transform of) the two- and four-legs vertex functions and at some 
special values of the external momenta characterized by some given scale s. Owing to 
our assumption of unbroken O(N) symmetry, we may choose to renormalize at zero 
momentum, while for r^ 4 ) we make the usual choice of the symmetric configuration. These 
special values of r^ 2 ^ and will be identified, respectively, with the physical squared 
mass at equilibrium and with the coupling constant renormalized at the scale s. 

Using the relations ([2.6)1 and definitions ([3.24)1 . ([3.25)1 we obtain the following form 
for the free mode functions 

The two-legs function is obtained by calculating ([3.30)1 in zero £ field. It is easy to check 
that V = is solution of eq ([3.30)1 for £ = if the mass renormalization eq. ([3.23)1 holds 
true. Using this fact and the general definition of vertex functions, eq. ([2.22)1 . we obtain 
the free two-points vertex function 

vf){x\y) = {U + m 2 )8 ij 8 A {x-y) 

We recall that this is a vertex function in the physical representation of the CTP formal- 
ism, with one a and one C insertions. 

Now we have to calculate the four-point function (one a and three <f) c ). From the 
variation of the second eq. ([3.30)1 with respect to £ at £ = we find the integral equation 



$ijkm(x) = A r ijkm 5 {4 \x) - ±A T ijrs I d A x' I {2 \x - x') d rskm (x') (3.32) 
for the function 
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where 

I {2 \*) = -^G% j { X )Gf{ X ) 

is the (suitably normalized) loop with one retarded and one correlation Green function. 
Eq. (|3.32|) becomes algebraic in Fourier space 

$ijkm(p) = A r ijfc m - |A I {2 \p) r ijrs $ rskm (p) (3.33) 

and is easily solved by 

%fc m (p;A ) = |#i(p;A ) [SikS jm + S jm 8 ih \ + ~# 2 (p; A ) 5ij$km (3.34) 

with gi and g 2 defined as 

9l(p;Ao) = i + iA A %» (p) ' * (p;A ° ) = T7§#T&M (X35) 

We recall that, as appropriate to a causal initial value problem, all Fourier transforms 
are analytic in the upper complex po^halfplane. Notice also that I^ 2 '(p) = I^ 2 \p,m 2 , A) 
depends on the UV cutoff A through the initial time correlation function G^\x), 

P\p,m 2 ,A) = -^j iki<A ^G^(p-k)G^(k) (3.36) 

so that, as long as A is finite, it is a function of both the time component po and of the 
Euclidean 3D length \p\, rather then just of the Lorentz invariant p 2 = p% — p 2 ■ 

We now apply the relation between V and T' in ()3.26j) and the general definition ()2.22|) 
and obtain the following form of the four-point vertex function 

r^ k j Xl \ X2 , X3 , X4 )= ! [n^]f^( P2 ,p S) p 4 ) c ^E.M----) 

J n=i W J (3.37) 

^jkm = _2 Ao njkm + $ijkm{P3 + Pi) + $ikjm{P2 + Pi) + $imkj{P2 + Pz) 

Notice the structure of as sum of three terms, for the s, t and u channels, with the 
same functional form given by ^ijkmip) ~ 2A /3. This structure is characteristic of the 
mean field approximation. 

In the present out-of-equilibrium context, we can define the symmetric point at the 
scale s as 

P3 + Pi = P3 + P2 = P2 + Pi = Qs = (0, q) , \q\ = s 

with an arbitrary direction q. Evaluating f^ 4 ) at this point at yields the renormalization 
condition 

Xr ijkm = f^ fcm | sympt At A = -2A + 2 gi (q s ; A ) + g 2 (q s ) A ) (3.38) 
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This result should be compared to the corresponding one in the static N = 1 case, that is 
eqs. (J3.7|) . (|3.11|) . We see that, in the current context with generic N and generic renor- 
malization scale, the effective coupling 3g is split into 2gi and #2, both scale-dependent. 

Before going any further in discussing renormalizability and RG invariance (indepen- 
dence on the scale s), it is convenient to calculate a "quasi-renormalized" form of the 
self-consistent eq. (|3.3U|) for T' analogous to the eq. (|3.1U|) found for the effective poten- 
tial. 

First of all let us define J[V) as the part of the integral term /[V] in eq. ()3.30|) that 
doesn't contain any divergent integral 

J{x) =I{x)- \G® + Gf V G<® + G<® V T Gf] {x, x) 

that is, in Fourier space 

I {2 \p)Vkm{p) (3.39) 

where 8^(p) stands for —i(po + ie)~ 1 5^(p), as appropriate in this out-of-equilibrium 
context. 

The second of eqs. ()3.3()|) can now be rewritten as 



Vij(p) = \$i jk m{p ; A ) [Lm{p) + J[V]km{p)} (3.40) 

Now, we can define 

A^r.(p) = P l3 {p) - \m 2 5^\p) - [\# ijkm (p) - lX 5 ik 6 jm ] £ km (j>) (3.41) 

and, by substituting in ()3.40|) . we obtain the HF self-consistency equation in the "quasi- 
renormalized" form that we were looking for, that is 



(3.42) 



which is to be compared with eq. ()3.10|) . 

We can now repeat the analysis of the previous subsection in this more general case. 
As we will see, conclusions will be rather similar. 

First of all, exactly as it happened in the static case, the renormalization of the mass 
parameter receives contributions form daisy and superdaisy tadpoles with the topologies 
shown in ()3.15|) . These contribution are actually not affected by momenta entering in 
the two-legs function and simply define the bare mass in terms of the renormalized one 
according to eq. (j3.23J) . From now on we deal with dressed propagators, i.e. propagators 
with the renormalized mass. 

We look then at the coupling renormalization in eq. ()3.38)) . It exhibits the same 
pathological behaviour of the static N — 1 case, eq. ()3.7|) . now with a iV— dependent 
"Landau obstruction". Most importantly, from the explicit form of the HF four-legs 
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function in eq. ()3.37|) we see that imposing a finite value to A at a certain chosen scale 
s fails to render finite the four-legs function with generic incoming momenta and, in 
particular, the coupling A' at any other scale s' ^ s. We can correct both shortcomings 
with the same strategy of the static N = 1 case. In fact, the diagrammatic analysis of 
the HF four-legs function shows that it is given by the resummation of "bubble" graphs 

+ V + s U m +\ ^'k^\ + s V m H 

' v - ' \ ' ^ - ' ^ - '\ (3.43) 

where we recall that dotted (solid) lines represent </>a (<fi c ) propagators and external legs 
coupled with the vertices (J2.19)) of the physical representation. This resummation misses 
all leading logarithmically divergent contributions coming from non bubble diagrams. In 
the usual in-out formalism their inclusion provides the standard 1LRG relation 

- ~ ~ + /3 A logA + ... 
A Ao 

where (3q = is the first coefficient of the beta function and the dots stand for the 
scheme-dependent finite parts. Now the HF approximation dictates a precise scheme 
where leading logarithms and finite parts must all be contained in the 1-loop bubble 
l( 2 '(q s , m 2 , A), implying the complete parametrization 

«A,^A/ n )| m = T -^-^ (3.44) 

which generalizes to the present non-static case with generic N the static N — 1 result 
in eq. (J3.17)) and replaces the pure HF parametrization of eq. (|3.38)1 . Notice that in the 
out-of-equilibrium formalism there could be a priori a complication with non-local finite 
parts, since the four-legs function even if evaluated at the symmetric point, need 
not be in principle completely symmetric in its four O(N) indices (recall that it has one 
0a leg and three <p c legs). However the symmetry in the three <p c legs is enough to fix 
the r tensor as the unique O(N)— invariant fourth-rank tensor, so that a unique coupling 
constant renormalization can occur, including finite parts, just like in the standard in-out 
formalism. 

We need now to determine how the new contributions which change the HF parametriza- 
tion of eq. (J3.38)) to the new one in eq. (J3.44)) modify the four-legs function away from 
the symmetric points of the external momenta. Consistently with the mean field nature 
of the approximation, r^ 4 -* must be the sum of three separated channels with the same 
functional form Ao/[Ao/ < - 2 - ) (p)] valid at the symmetric point. Then the request of finiteness 
plus agreement with the perturbative one-loop result uniquely fixes the following specific 
form 

ikjm ijkm 

(3.45) 

where the tensor I reads 

1 



-ijkm 



3(iV + 8) L 



(N + 4)5ij8 km + 2(5 ik 5 jm + 5 im S jk ) 
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and X r {p) is the running coupling constant, connected to coupling A = X r (q s ) at the scale 
[i by the relation 

_L_4 + ^ [ /« W _/« M ,3.46) 

that involves by construction only a subtracted integral, so that one may now take the 
limit A — > oo. Of course, this defines a meaningful X r (p) only for p timelike or p spacelike 
but such that — p 2 is smaller than the Landau pole. 

At any rate, eq. ()3.45|) shows explicitly that requiring a finite coupling strength at a 
scale s makes the whole four-legs function finite. Furthermore, the expression ()3.46j) is 
by construction RG invariant. Therefore it seems sensible to require such a form for the 
coupling and modify the HF resummation by including the missing contributions that 
appear in the CTP diagrams with the topology in eq. (j3.17|) . 

We still have to consider the higher order vertex functions. These are encoded by 
construction in the functional A.7^ defined in (I3.41j) . From the HF self-consistency 
equation ()3.42|) for AjF^. we can see that an explicit dependence on the logarithm of the 
cut-off is still present even after the coupling constant renormalization. The reason is 
in the presence of the effective vertex dij km (p) in eq. ()3.4()|) . which corresponds to CTP 
diagrams of the same type of those drawn in (j3.19|) . Exactly in the same way of the N — 1 
static case [see eqs. f!3.18|) and (j3.19p ] we have to include other contributions containing the 
correct leading logarithms that come from the 4>a — 4>c diagrams corresponding to those 
in eq. ()3.20|) and (j3.21|) . In order to obtain renormalizability and maintain the single 
channel structure of the effective vertex (which is a direct consequence of the mean-field 
approximation), these contribution have to be taken momentum dependent in such a way 
to determine a Leading Log structure for the effective vertex in place of $ij km (p). Then 
finiteness and agreement with the lowest pertubative order require the resulting structure 
to be simply X r (p) r ijkm . 

In conclusion, we are naturally led to consider an extended resummation of dia- 
grams that will satisfy the requests of renormalizability and RG invariance, by modifying 
eq. ()3.42|) in the following way 

r l3 {p) = Ai* (p) + > 2 M (4) (p) + \K(p) tmmZUp) 

A.^0) = \K{p) r ijkm J[V] km (p) (3.47) 

%{p) = 2APlj(p) + lK{p)njkmikm(p) 

The renormalized equations of motion are then 

[(□ + m 2 )(% + Vij(x) - lEij(x)] ^(x) = 

[(□ + m 2 )5ij + Vij{x)] u kaJ (x) = (3_ 48 ) 

%(P) = |A r (p) T ijkm ikm(p) + hm(p) ~ I < E ) S ( ' 4) (p)5 km + J (2) (p) V km (p) 

where ^ij(x) is the inverse Fourier transform of H^(p) = §A r (p) {jij km — iij km ) Ckm{p) and 



23 



we recall that £(x) = (pi(x)(pj(x) and hmip) ls the Fourier transform of 



I[V} lj (x) = Re 



d 3 p 



p\<A 



(27r) 3 Upb ' i ^ x ^ u P b >i\ x ) 



which is a functional of the mean field V through the evolution equation for the mode 
functions. 

The absence of any divergent cut-off dependence in the modified HF approximation 
defined by eqs. (|3.47j) or by the equations of motion (J3.48j) is manifest, since they involve 
only the finite running coupling constant and the fully subtracted (Fourier transform of) 
the equal point correlation hj{x). Comparing the expression for Vij(p) with the original 
one in eq. (j3.30j) . that reads in Fourier space, 



Vijip) 



1 \ 

2 'ijkn 



lknip)+hn{p)-I i E ) ^\p)5 kn 



we see that our improved renormalization simply amounts to perform the logarithmic 
subtraction on /[V]y while promoting the bare coupling Ao to the running one X r (p). As 
a consequence, the p-dependence of A r (p) and I^ 2 \p) imply a space-time nonlocality of 
the modified HF self-consistency equation for the mean field V(x), which now depends 
on the whole history of the background (f) and of the mode functions Uk a from time zero 
to time t. A similar space-time nonlocality appears in the equation of motion for the 
background field 0. Causality in this nonlocal evolution is guaranteed by analyticity in 
the upper p — halfplane of I^(p). Let us observe, however, that the third equation in 
eqs. ()3.48j) does not provide an explicit expression of the mean field V in terms of the 
background and of the mode functions. To obtain such an expression, we must solve for 
V, paying the price of losing manifest finiteness. In practice, this reduces to the simple 
matrix inversion already performed on eq. ()3.32|) and yields 

%(P) = -l$ijkm(p;-K(p)) [ikm(p) + hm(p) ~ 1$ 5® (p)5 km \ 

where the fiijUmip \ u) is the function defined in eq. ()3.34|) . In coordinate space this becomes 



(i)' 



Vij( x ) = 3 J d x |2 7i(x - x ) 4>i(x')4>j(x') + hj(x') - SijI E 
+ l2 (x - x') U 2 (x') + I kk (x') - iv4 1} 



(3.49) 



where ji(x) and 72(2;) are the inverse Fourier transforms of the functions — <7i(p; — A r (p)) 
and —g2(p', —Kip)) 7 respectively, which can be read from eq. (|3.35|) and are indeed ana- 
lytic in the upper po— halfplane. For instance, we have explicitly 



7i {x) 



d A p 



A e~ ip -* 



(2vr) 4 1 + 1A[(JV + l)K 2 )(p) - (N + l)lW(q s )] 



with an even more involved expression for 7 2 (x). Clearly 7 12 depend logarithmically 
on the cut-off A through an incompletely subtracted I^(p), but, by construction, this 
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logarithms cancels out with those coming from the incomplete subtraction hj{x) — Sijl^ 
on the equal point correlation of the mode functions. The memory functions 71 (x) and 
72 (x), as well as the inverse Fourier transform of \ r (p) which enters in the evolution 
equation of the background field [see eq. (|3.48j) ]. all share the same behaviour for large 
timelike values of x, namely a (rr 2 )~ 1//2 decay modulated by oscillating factors. 
When N = 1 things simplify considerably and we find for the mean field 



V(x) 



J d 4 x' j(x — a/) j 



</>V) + 



d 3 k 
(2^)3 



where 



1 (x) = §7! (re) + 572O1O 



d 4 p 



\u k (x 



X e~ ip - 



2u(k) 



3,u ' 3 ' ZK ' 2 ' (2vr) 4 l + A[/( 2 )(p)-|/( 2 )(g, 
while mode functions and background evolve according to 



[□ + m 2 + V(x)] u k (x) = 
□ + m 2 + V(x) - § / d 4 x' 7 r (x - x')(j) 2 {x') <f>(x) = 



where 



J r (x) 



d A p 



Xe 



■ip-x 



(27r) 4 l + |A[/( 2 )(j9)-/( 2 )(g s 

is the ultraviolet finite inverse Fourier transform of the running coupling constant X r (p). 

On the opposite end, when iV — > 00, one can verify that the standard local evolution 
equations [H] are recovered from eqs. ([3.48)1 . To this end, it is convenient to restrict 
first to the special case of a background field which maintains a fixed direction, that is 
4>i(x) = <j)(x)vi, with v some fixed N— dimensional unit vector. In fact in this case we 
can reduce eqs. ()3.48|) to an index-free form by introducing longitudinal and transversal 
projectors 



L,ij 



ViVj 



T,ij 



'ij 



L,ij 



Then we can set 



Uka,i( X ) = u kL(x)PL,ia + U kT (x)P T ,ia 

Vij(x) = V L (x)P L>ij + V T (x)P T4j 

which are compatible with initial condition in eq. (j3.31|) if we assume that the longitudinal 
and transverse parts of the mode functions and their derivatives are equal at initial time. 
By substituting into the equations of motion and projecting one obtains 








[□ + m 2 + V T (x)}u kT (x) 
[D + m 2 + V L (x)}u kL (x) 
[□ + m 2 + V L {x) - \<b(x)\<j>(x) = 

MP) = IHp) + lK{p) [(N - l)Jr(p) + SJ L (p) 

V T (p) = i$(p) + |A r (p) \(N + 1) J T (p) + J L (p) 



(3.50) 
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where $(x) is the inverse Fourier transform of $(p) = A r (p) £(p) [recall that £(x) = 4> 2 {x)} 
and the meaning of Jl and Jt is obvious [see eq. ()3.39|) ]. Again, the formulation in 
momentum space for the mean field V, while convenient to avoid lengthy convolution 
integrals and to allow for an easy check of finiteness, does not make causality manifest. 
The opposite situation, with finiteness hidden and causality manifest would follow by 
solving explicitly for V and then reverting to coordinate space as done above for the 
general case [see eq. ()3.49|) ] or for the N — 1 case. 

Now let us rescale the coupling and the background field as prescribed by the standard 
large N procedure 

A -> X/N <f>(x) -> y/N<p{x) 

By substituting in ()3.50|) and taking the limit on N we obtain (notice that the longitudinal 
mode functions decouple) 



[□ + m 2 + V T {x)]u kT {x) 
[□ + m 2 + V T (x)](f>(x) = 



V T (x) = ±A 4>\x) + 



= 

= 

d 3 k 
(2tt)3 



1 



u k (x 



T (X, 



2co(k)l 

which are the usual, local and renormalized large N evolution equations. We remark that 
A is the renormalized coupling at the scale s while the nonlocal terms (with the running 
coupling A(p)) have all disappeared by explicit cancellation. This is indeed what was to be 
expected, on the basis of the same diagrammatic analysis that lead to our improvement of 
the HF approximation, since the nonrenormalizable terms in the standard HF approach, 
as well as the contribution we have added to cure the problem, all come from diagrams 
that are at least 1/N suppressed. 

We conclude this lengthy subsection with an explicit check of the Renormalization- 
Group invariance of our improved HF approximation, by verifying that indeed all vertex 
functions obtained by eqs. (|3.47|) are solutions of Callan-Symanzik equation. These equa- 
tion states the RG invariance of any observable O which is a function of coordinates 
(or momenta Pi), of the scale s, the coupling A and the parameter a = m 2 /s 2 , that is 



s-^+(3(X,a)-^ + f3^(X,a) ° 



ds 



da 



0({x t }\s,X,a) = 



(3-51) 



with {3 and (3 $2 functions to be determined. Notice the absence of the term with the 
anomalous dimension due to the lack of field renormalization in HF approximation. By 
applying eq. (J3.51|) to the two-legs function one obtains (3^2 = —2a. Then by applying to 
four-legs function one obtains 



/5(A,a) = -|(iV + 8)A 2 



lim 

A— too 



m 



A) 



s 2 m 2 we have (3 ~ A 2 as expected. Now, is easy to check that all free 



Then for 

propagators with dressed mass satisfy the CS equation since they are functions of m 2 
alone. Then, since the generic ra-legs function is a functional of the propagators and of 
the RG-invariant running coupling \ r (p) through the four-legs function, one immediately 
concludes that it satisfies eq. (|3.51|) . 
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3.3 Other initial states and nonzero temperatures 

As stated in the previous subsection, with generic initial conditions on the background 
field eqs. (I3.48J) already describe an out-of-equilibrium problem, in spite of the choice in 
eq. (|3.31|) of equilibrium initial conditions on the quantum fluctuations. It is nonetheless 
sensible to ask whether and how we can choose different initial conditions for the mode 
functions without spoiling the properties of renormalizability and RG invariance. 

What might happen can be shown by the following example. Let us consider, in the 
simple N = 1 case, an initial state of the same form of the HF vacuum but with a different 
mass M, which could be for instance the solution of the gap equation ()3.3|) in the case 
of an uniform background. This is a frequent choice in dealing with non-equilibrium 
problems and has a precise physical meaning, since it corresponds to the minimization 
of the quantum fluctuation part of the HF energy at fixed uniform background. At any 
rate, as soon as we assume as initial conditions 

Si 



with M m, then the free mode functions have more than one frequency component 



)y/2n(k) I 



Lj{k) 



e -iu{k)t _|_ 



f(fc). 



e -iw(fe)t 



and the free correlation function is no longer translationally invariant in time. Then the 
integral term in eqs. (|3.48j) [recall that here N = 1] 

in spite of the subtractions still contains the superficially divergent contribution 

f d 3 k m 2 -M 2 . K , 

that indeed diverges with the cut-off when t = 0. These initial time singularities have 
first been discussed in [35. and removed by a Bogoliubov transformation on the initial 
state, which in practice amounts to a redefinition of the initial kernel in such a way that 
the leading terms of an high-momentum expansion are the same as at equilibrium 

^ 1 m 2 

y(k) ~ — = + — . + . . . 
2 v / fc 2 4(fc 2 ) 3 / 2 

In fact, one may verify that the initial singularity (and any other divergence as well) are 
absent for any choice of kernel having the above large k expansion. A simple interpretation 
is that the renormalization procedure ensures finiteness for any initial Gaussian state 
belonging to the same Fock space of the HF vacuum. From this simple example we 
can extrapolate the generic condition on the short-distance behaviour of the initial state 
kernel 

1 m 2 

G(x, y) ^ -r-2] r + ^-o lo s \ x - v + • • • 

4ir z \x — y\ 8tt 2 
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that ensures the cancellation of all divergent terms are guaranteed by mass and coupling 
constant renormalization. This conclusions can be immediately extended to the generic 
N > 1 case. 

Before concluding let us say some words about the case when the initial state has a 
nonzero temperature T. The formalism for pure state dynamics introduced in section |2~21 
can be easily generalized to statistical mixtures defined by Gaussian density matrices. 



p[(pi,(p 2 ) =A/"exp ji (p\<p! - Lp 2 ) - (A<pi| [\Q 1 +iS] |A^i) 

- (Ap 2 | [ig- 1 - is] \a V2 ) + \ (A^| g-\ ^g-\ \a V2 ) } 



(3.53) 



Where Atpi = tpi — <fi. The parameters are the the background field <fi and its momentum 
p, the symmetric kernels g, S, 7 and the antisymmetric a (with C = 7 + 4i£ 1//2 a g 1 ^ 2 ). 
The normalizability of the state (Tr[p] = 1) requires g~ x l 2 (1 — 7) g^ 1 ^ 2 to be positive. 
In CTP formalism the generalization proceeds simply by substituting ^[<p + ]^'[(p_] with 
p[<p+, <p-\ in the path integral of eq. (|2.10p . Then HF equations are still as before while 
more general initial conditions for the mode functions are allowed. Let us consider an 
initial state with translationally invariant kernels and, for simplicity, let us put N — 1, 
S = a = 0, then we have 

Uk (x,o) = [g(k)]^ 2 [i-m. 



-1/2 ikx 



ii k (x,0) 



^(fc)- 1 [l-f^)] 1 / 2 ! u k (x,0) 



For vanishing initial background field and momentum, standard equilibrium solutions are 
given by 



g(k) 



tanh 



2u T (k) T 
with mtp defined by the gap equation 



LU T (k) 



k 2 



m 



T 1 



Ik 



cosh 



L)T{k) 



T 



(3.54) 



rrin 



m 



./W + (m | _ m 2 )/( 2) (0)+ l A J 



d 3 k 



1 



coth 



u T (k) 



(2tt) 3 2u T (k) 2T 

By expanding ujT{k) around uio{k) = u(k) one can easily check the cancellation of diver- 
gent terms. We can then consider the out-of-equilibrium problem with the same initial 
conditions on the mode functions but arbitrary initial values for the background field and 
its velocity. The free mode functions read 

1 



a 



.r 



coth ^ 

2lu t 2T L 



1 + 



to 



-iujt 



LU 



Axk 



Calculating the free correlation function we easily can check that the mean field expression 
in the last of eqs. ()3.48j) is free of divergences, except for the one at t — 0, which is given 
by eq ()3.52j) with M = tut- As above, this can be cured by a Bogoliubov transformation, 
which in practice amounts to changing the equilibrium solution in eq. ()3.54|) by allowing 
a suitable ultraviolet fc-dependence in the temperature. 

One last comment regards RG invariance. We can see that the proof of the subsec- 
tion 13.21 still holds as long as the free correlation function depends on the renormalized 
squared mass alone. 
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